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Abstract 

Q_i' We use on-shell recursion relations to determine the one-loop QCD scattering amplitudes with 

a massless external quark pair and an arbitrary number (re — 2) of positive-helicity gluons. These 



amplitudes are the last of the unknown infrared- and ultraviolet-finite loop amplitudes of QCD. 
The recursion relations are similar to ones applied at tree level, but contain new non-trivial features 
corresponding to poles present for complex momentum arguments but absent for real momenta. We 
present the relations and the compact solutions to them, valid for all n. We also present compact 
forms for the previously-computed one-loop n-gluon amplitudes with a single negative helicity and 
the rest positive helicity. 
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I. INTRODUCTION 



The computation of new one-loop amplitudes in perturbative gauge theories, and in 
QCD in particular, will be a prerequisite for theoretical studies related to the experimental 
program at CERN's upcoming Large Hadron Collider (LHC). The discovery and study of 
new physics beyond the standard SU(3) x SU(2) x U(l) model of particle interactions will 
depend on our ability to calculate higher-order corrections to a wide variety of processes in its 
component gauge theories. Computations of tree-level scattering amplitudes are a first but 
insufficient step. The size and scale-variation of the strong coupling constant imply that a 
basic quantitative understanding must also include the one-loop amplitudes which enter into 
next-to-leading order corrections to cross sections. Such corrections are also required to build 
a theoretical base for a program of precision measurements at hadron colliders [1] . Precision 
measurements at the SLAC Linear Collider (SLC) and CERN's Large Electron Positron 
(LEP) collider have proven the power of such a program in advancing our understanding of 
short-distance physics. 

Recent years have seen important progress in this theoretical program; yet a wide and 
seemingly hostile province still severs us from our goal, encouraging us to seek additional 
tools for performing the required calculations. The last year in particular has seen new 
progress in the computation of tree-level [2, 3, 4, 5, 6, 7, 8, 9, 10] and one-loop [11, 12, 13, 
14, 15, 16] gauge-theory amplitudes stimulated by Witten's proposal [17] of a weak-weak 
coupling duality between Af = 4 supersymmetric gauge theory and the topological open- 
string B model in twistor space, generalizing Nair's earlier description [18] of the simplest 
gauge-theory amplitudes. Further investigations along these lines have revealed new aspects 
of the underlying twistor structure of gauge theory [19, 20, 21, 22, 23, 24, 25]. For a recent 
review, see ref. [26]. 

Recursion relations, originally introduced by Berends and Giele [27], have long been recog- 
nized as an efficient and elegant method for computing tree-level amplitudes. Other related 
approaches [28], as well as various computer-driven approaches such as MADGRAPH [29], 
have also been employed. Recently, Britto, Cachazo and Feng wrote down [6] recursion 
relations, employing only on-shell amplitudes (at complex values of the external momenta). 
These relations were stimulated by the compact forms of seven- and higher-point tree am- 
plitudes [5, 13, 15] that emerged from infrared consistency equations [30]. They were shown 
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to yield compact expressions for next-to-maximally-helicity-violating (NMHV) tree ampli- 
tudes [8]. The same authors and Witten gave a simple and elegant proof [7] of the relation 
using special complex continuations of the external momenta. 

The proof, which we review in section III, is actually quite general, and applies to any 
rational function of the external spinors satisfying certain scaling and factorization prop- 
erties. Indeed, it has since been applied to amplitudes with massive particles [10], and in 
gravity as well [31]. The generality of the proof suggested that it should be useful for finding 
on-shell recursion relations at one loop. We previously wrote down [32] such relations for 
the (infrared and ultraviolet) finite n-gluon amplitudes in QCD, An\l ± , 2 + , 3 + , . . . , n + ), for 
which all gluons (or all but one) have the same helicity. Unlike the situation for massless 
tree amplitudes, in an application to general loop amplitudes in a non-supersymmetric the- 
ory, factorization in complex momenta is qualitatively different from that in real momenta. 
Accordingly, we had to address new issues, in particular the appearance of double poles in 
the complex analytic continuation. 

In this paper, we examine another application of such on-shell relations, to one-loop 
n-point amplitudes with one pair of massless external quarks and (n — 2) positive-helicity 
gluons, An\l$ , 2+, 3 + , . . . , n + ). This set of helicity amplitudes, together with the above 
n-gluon amplitudes and their partners under parity, are zero identically at tree level due to 
supersymmetry Ward identities [33]. This is because massless quarks differ from gluinos only 
in color manipulations which are essentially trivial at tree level. At one loop, the difference 
in color factors between quarks and gluinos permit the amplitudes to be non-vanishing. 
However, any infrared and ultraviolet divergences would have to be proportional to the 
corresponding tree amplitude, which vanishes. Hence this set of helicity amplitudes is finite. 
Because all the "zeroes" have been filled in at one loop, none of the corresponding two- 
loop helicity amplitudes can be finite. For example, the first two-loop four-gluon scattering 
amplitude to be computed [34], A\ (1 + , 2 + , 3 + , 4 + ), has infrared divergences similar to a 
typical one-loop amplitude. Thus the amplitudes we compute in this paper represent the 
last finite loop amplitudes to be computed in massless QCD. 

We calculated the five-point amplitude, An\lg, 2+, 3 + , 4 + , 5 + ) (together with all the 
other helicity configurations) long ago [35], but no other results in this class of amplitudes 
exist in pure QCD. On the other hand, related QED and mixed QED/QCD amplitudes, 
containing a massless external e + e~ pair and arbitrarily many positive-helicity photons or 



3 



gluons, were computed by Mahlon [36, 37]. 

In the recursive approach to the finite quark-gluon amplitudes that we take, we must 
treat a new type of pole in addition to the double poles encountered in the pure-gluon 
case [32]. These are poles where the collinear limit in real momenta is not singular, yet 
nonetheless a single pole arises for complex momenta. We shall call these "unreal" poles. An 
understanding of these unreal poles is essential for constructing correct recursion relations. 
In this paper we will not provide a first principles derivation of the unreal poles, but will 
instead take a pragmatic approach. First we determine the unreal pole contributions to the 
recursion relation by reconstructing the known five-point amplitudes [35]. Then we use our 
determination of the unreal poles at five points to write down a pair of recursion relations 
for an arbitrary number of external legs. We find compact solutions to the two relations, 
valid for all n. For a subset of the amplitudes, we also find a set of recursion relations in 
which unreal poles are absent. 

To confirm these relations we perform non-trivial consistency checks of the factorization 
properties of the solutions. As an additional cross check, we use our QCD amplitudes to com- 
pute mixed QED/QCD and pure QED amplitudes, by carrying out sums over appropriate 
permutations of the gluon momenta. The permutation sums are designed to remove non- 
abelian self-couplings, thus allowing a quark-anti-quark pair to mimic an electron-positron 
pair. We then compare the permutation sums with Mahlon's earlier results for the same 
amplitudes [36, 37]. 

As a side benefit of our analysis we also find a compact form for the previously com- 
puted [37] one-loop n-gluon amplitudes with a single negative helicity and the rest positive. 
This form is obtained from the quark amplitudes by taking the momenta of the quark and 
anti-quark to be collinear, and making use of the known [38, 39, 40] factorization properties. 

On-shell recursion relations may provide a technique, complementary to the unitarity- 
based method of Dunbar and the authors [34, 38, 41] (and its more recent refinements [14]), 
for the calculation of one-loop QCD amplitudes. The unitarity-based method applies most 
easily to terms in the amplitudes that have discontinuities. Computation of these terms 
requires only knowledge of tree amplitudes evaluated in four dimensions. The unitarity- 
based method can also be applied to computation of rational terms, by computing the cuts 
at higher order in e [34, 41]; but doing so requires knowledge of tree amplitudes with two legs 
in D = 4 — 2e dimensions. These are equivalent to amplitudes with massive scalars [41, 42], 
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to which on-shell recursion relations have also been applied recently [10]. 

This paper is organized as follows. In section II, we introduce the notation used in the 
remainder of the paper, and describe the color organization of one-loop QCD amplitudes 
with an external qq pair, in terms of color-stripped objects called primitive amplitudes. In 
section III, we review the general construction of recursion relations, as well as summarize 
the known factorization properties of one-loop amplitudes, which dictate the structure of 
the relations. The recursion relations for the quark amplitudes are built using a known 
set of amplitudes which we present in section IV. In section V, we construct our relations 
for the finite quark amplitudes, solve them, and discuss the factorization properties of the 
solutions. We provide additional cross checks on the solutions in section VI. 

II. NOTATION 

We will organize the calculation of one-loop two-quark (n — 2)-gluon amplitudes in terms 
of primitive amplitudes [35] . Each color-ordered amplitude in a trace-based color decompo- 
sition [35, 43, 44, 45, 46, 47] is built out of several primitive amplitudes. 

For tree-level amplitudes with a quark pair in the fundamental representation and (n — 2) 
external gluons, the color decomposition is [44, 45, 46], 

A^(l q ,2 q ,3,...,n) = g n ' 2 £ (TM«> . . . TM»>)* A£\l„ 2 q ; <j(3), . . . , a(n)) , (2.1) 

creSn-2 

where S n -2 is the permutation group on n — 2 elements, j hj denotes the j-th momentum 
kj and helicity hj, and the superscript "(0)" signifies that these are leading-order tree am- 
plitudes. The T a are fundamental representation SU(iV c ) color matrices normalized so that 
Tr (T a T b ) = 5 ab . The color-ordered amplitude An^ is related to tree- level all-gluon ampli- 
tudes by supersymmetry Ward identities [33]. Because the color indices have been stripped 
off from the partial amplitudes, there is no need to distinguish a quark leg q from an anti- 
quark leg q; charge conjugation relates the two choices. 

At one loop an additional trace of gluon color matrices T ai may survive. The general 
color decomposition for fundamental representation quarks is [35], 

A(l 5 ,2„3,...,n) = tJ^E E Gv^(a(3,...,n))A n . J (l q ,2 q ;a(3,...,n)), 

{ ' 3=1 aeS n - 2 /S n;j 

(2.2) 
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where we have extracted a loop factor of l/(Air) 2 , and the color structures Grj^ are defined 
by, 

N c (T aa . . . T a ") J 1 , 

, (2.3) 
Tr(T a3 . . . T a ^ +1 ) (T a ^+ 2 . . . T a ") J 1 , j = 3, . . . , n - 2, 
Tr(T a3 2~ ,a ") 

Here S^j = Z,-_i is the subgroup of SVj-2 that leaves Gr^ invariant. When the permutation 
a - acts on a list of indices, it is to be applied to each index separately: <r(3, ...,n) = 
<r(3), . . . , cr(n), etc. We refer to A n; i as the leading-color partial amplitude, and to the 
A n ;j>i as subleading-color, because for large N c , A n -i alone gives the leading contribution 
to the color-summed correction to the cross section, obtained by interfering An^ with A n . 
The explicit N c in the definition of the leading-color structure Gr^ — which is otherwise 
identical to the tree color structure — ensures that A n;1 is 0(1) for large N c . (For super- 
Yang-Mills theory, where the fermions are adjoint-representation gluinos, one should use the 
same color decomposition as for gluons [47]. However, the particular helicity amplitudes 
considered in this paper vanish for super- Yang-Mills theory.) 

We describe the amplitudes using the spinor helicity formalism [46, 48]. In this formalism 
amplitudes are expressed in terms of spinor inner-products, 

(jl) = (r\l + }=u4k 3 )u + (k l ), [jl] = <j + |/-> = «+(*,-)«-(*,) , (2.4) 

where u±(k) is a massless Weyl spinor with momentum k and plus or minus chirality. Our 
convention is that all legs are outgoing. The notation used here follows the standard QCD 
literature, with [ij] = sign(kfk'j) (j i)* so that, 

(ij)[ji] = 2ki- kj = Sij . (2.5) 

(Note that the square bracket [ij] differs by an overall sign compared to the notation com- 
monly used in twistor-space studies [17].) 

We denote the sums of cyclicly-consecutive external momenta by 

K.. 3 = fcf + k f +1 + .-. + k^ + fcj, (2.6) 



Gr^(3,...,n) = 

Grg } (3;4,...,n) = 

Grg ) (3,...,j + l;j + 2,...,n) = 

Grg?-i(3,...,n) = 
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(a) (b) 

FIG. 1: In diagram (a) the fermion line (following the arrow) turns left on entering the loop and 
is an L type primitive amplitude. In diagram (b) the fermion line turns right and is an R type. 

where all indices are mod n for an n-parton amplitude. The invariant mass of this vector is 
Sj...j = Kf...j- Special cases include the two- and three-particle invariant masses, which are 
denoted by 

s ij = j = i^i kj) = ■ kj, Sijk = {ki + kj + kk) . (2-7) 

In color-ordered amplitudes, only invariants with cyclicly-consecutive arguments need ap- 
pear, e.g. Si t i + i and Sj,j+i,j+2- We also write, for the sum of massless momenta belonging to 
a set A, 

K^Y. K ■ (2-8) 

Spinor strings, such as 

(i + \$ A \j + ) = Y,[^](*3) , (i-\$A\r) = Y,(ia)[aj] , (2.9) 

aeA a£A 

and 

(r\(a + b)\j-) = (ta)[aj] + (tb)[bj} , 

(r\(a + b)(c + d)\j + ) = (ia)(a + \(c + d)\j + ) + (ib)(b + \(c + d)\j + ) , (2.10) 

will also make appearances. 

As noted above, we can write the one-loop color-ordered amplitudes in the fermionic case 
more compactly in terms of primitive amplitudes. Each primitive amplitude corresponds to 
the set of all color-ordered diagrams with specified internal states, and a specified orientation 
of the fermion line along the loop [35]. 

First consider the case depicted in fig. 1, where the fermion line entering the diagram 
is part of the loop. Upon entering, the fermion line can turn either "left" and circulate 
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(a) (b) 



FIG. 2: In diagram (a) the external fermion line passes to the "left" of the loop, following the 
fermion arrow, and is designated an L type. In (b) it passes to the "right" and is an R type. A 
gluon, fermion or scalar can circulate in the loop. The same decomposition also holds even if we 
emit additional gluons off the external fermion lines. 

clockwise, or turn "right" and circulate counter-clockwise around the loop. The two orien- 
tations correspond to separate primitive amplitudes which we designate as "L" and "R". 
They would carry the same color factor were the fermions in the adjoint representation, but 
carry different ones when the fermions are in the fundamental representation. This division 
is therefore gauge invariant. The second case, shown in fig. 2, is where the external fermion 
line does not enter the loop. This case may also be divided into "left" and "right" pieces, 
which we again label by "L" and "i?" . Here the division is based on whether the loop passes 
to the left or right of the fermion line. 

The primitive amplitudes for qqgg • ■ ■ g are, 

^ [J1 (l/,3,4,...,i, 2/ ,i + l,...,n), i q 

A«' [J] (l f ,3,±,...,j,2 f ,j + l,...,n), ' ' 2 ' ' 

where J = \ and J = denote the contributions with a closed fermion loop and closed 
complex scalar loop, respectively. The second fermion has been placed in the j th position, 
but assigned the label 2. Since the primitive amplitudes can be used to build amplitudes 
with any color representation for the fermions, instead of labelling the fermionic legs by q 
and q we label them by / to denote a generic fermion in any color representation. The 
normalization is such that two helicity states (Weyl fermions or complex scalars) circulate 
in the loop. Diagrams without closed fermion or scalar loops are assigned to J = 1; they 
may or may not contain a closed gluon loop, as the two types of diagrams mix under gauge 
transformations. For notational simplicity, we shall suppress the superscript "[1]", 

A^A^, A2 = A*W. (2.12) 

The primitive amplitudes (2.11) are not all independent. The set of diagrams where the 



incoming leg 1 turns left is related (up to a sign) to a corresponding set where it turns right. 
This relation is via a reflection which flips over each diagram, reversing the cyclic ordering: 

A*M(l f , 3, 4, . . . , 2 f , . . . , n - 1, n) = {-l) n A L J J \l f , n, n - 1, . . . , 2 f , . . . , 4, 3). (2.13) 

In addition, the super- Yang-Mills partial amplitudes for two gluinos and n— 2 gluons A^ USY = 
A^ SY are given by the sum (with all cyclic orderings identical) 

.SUSY _ /i SUSY — aL , aR , A L,[l/2] , jfl,[l/2] (9 u \ 

/i n — /i n;1 — A n -f- A n -f- A n -f- /i n , \ Z - Lq: ) 

because the "left" and "right" diagrams have the same group-theory weight for an adjoint- 
representation fermion. Due to supersymmetric cancellations between "left" and "right" 
primitive amplitudes, A| USY is always simpler than either or A„. Equation (2.14) allows 
one to obtain one of the four terms on the right with no effort, given A^ USY . In the case of 
the finite amplitudes we are considering in this paper, A^ USY vanishes. We will choose to 
eliminate A^, and compute A%. 

Finally, the following fermion-loop contributions vanish, 

A*^(l f: 2 f ,3A,...,n) = ^ 1 / 2 l(l / ,n,...,4,3,2 / ) = 0, 
AW](i / ,3,2 / ,4,...,n) = ^ 1 / 2 l(l / ,n,...,4,2 / ,3) = 0, (2.15) 

and similarly for the scalar-loop contributions. The restriction to "left" or "right" diagrams 
combines with the ordering of the external legs to leave only tadpole and massless external 
bubble diagrams behind; but these are zero in dimensional regularization. 

For scalars and fermions in the fundamental representation circulating in the loop, the 
leading-color contribution to eq. (2.2), A n;1 , is given in terms of primitive amplitudes by, 

A n;1 (l,-,2 ? ;3,...,n) = A L n {l f , 2 f , 3, . . . , n) - ^^(l /; 2 f , 3, . . . , n) 

+^^' [1/2] (V, 2/, 3, . . . , „) + ^4£[°] (1,, 2/ , 3, . . . , n) . (2.16) 

For QCD the number of scalars vanishes, n s = 0, while Uf is the number of light quark 
flavors. (Our fundamental representation scalars are normalized so that they carry the same 
number of states as Dirac fermions, 4iV c .) 

As explained in ref. [35], the subleading-color partial amplitudes A n;j>1 appearing in 
eq. (2.2) may be expressed as a permutation sum over primitive amplitudes, 

A n .j(lq, 2 q ; 3, . . . , j + 1; j + 2, j + 3, . . . ,n) 
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(-I)'" 1 J2 ^ 1] Kl/,2/,3,...,n)) 

aeCOP{a}{/3} L 

-2j-A*M (a(l f , 2/ , 3, . . . , n)) - («,(!,, 2 /? 3, . . . , n)) 



, (2.17) 



where a< G {a} = {j + 1, J, ... ,4, 3}, fa e {(3} = {1,2, j + 2, j + 3, . . . , n - l,n}, and 
C0P{a}{/3} is the set of all permutations of {1,2,..., n} with leg 1 held fixed that preserve 
the cyclic ordering of the ccj within {a} and of the fa within {/?}, while allowing for all 
possible relative orderings of the a* with respect to the fa. For example if {a} = {4, 3} and 
{/?} = {1,2,5} (the case required for A 5 . 3 ), then COP{a}{{3} contains the twelve elements 

(1, 2, 5, 4, 3), (1, 2, 4, 5, 3), (1, 4, 2, 5, 3), (1, 2, 4, 3, 5), 

(1,4,3,2,5), (1,4,2,3,5), (1,2,5,3,4), (1,2,3,5,4), (2.18) 

(1, 3, 2, 5, 4), (1, 2, 3, 4, 5), (1, 3, 4, 2, 5), (1, 3, 2, 4, 5). 

Thus, all partial amplitudes appearing in eq. (2.2) are expressible as sums over primitive 
amplitudes and it is sufficient to compute the primitive amplitudes in order to fully specify 
the complete color-dressed amplitudes. 

For the case at hand, where all external gluons carry the same helicity, a supersymmetry 
identity [33] may be used to prove that the fermion loop and scalar loop are the same up to 
a sign, 

^(17,2+,...,^,...,^) = -^1(^,2+,...,^,...,^) 

= .l^l;^' .// //•). (2.19) 

Thus, by computing the closed scalar-loop primitive amplitude we obtain also the closed 
fermion-loop primitive amplitude. 

We shall find it convenient to compute the combinations 

K, and A L n - s = A L n -A s n , (2.20) 

instead of A s n and A%- 

Furthermore, combining reflection symmetry and supersymmetry leads to, 

A L n (l-, >■ , f n+) = (-1)^(17, n + , ...,£+..., 2+) 

= (-1)" k USY (i/-, n+, .... Jf 2-)- A L n (iy, „■ , f 2 + ) 
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- 4^ (!-,„+, ....jf 2")- ^^(i", n +, . . . , £ 

(-lrf-^V" n + ,...,j> + ,...,2 + ) 

+ 4^(1/, n+ ....J/. ....2-)- 4^(1/, n+ ..... /f ... 
(-lrf-^V- n + ,...,j> + ,...,2+) 




2 + ) 



(2.21) 



+ A* (17, //• ./ r 2+) + (-1)^(17, 2 + , . . . , j+ 




so that 



^-(17, 2- , f n + ) = (-I) n+1 ^-(17, n+, . . . , j+ . . . , 2+), (2.22) 



which we can use to obtain A%~ s for j > \{n + l)/2] (\x~\ is the smallest integer greater 
than or equal to x). Note also that, as discussed above, A s n vanishes if j — n or j — n — 1. 

In summary, for amplitudes with a single quark pair with identical helicity gluon legs, 
there are two independent classes of primitive amplitudes that need to be computed, 



A L n ~ s {lJ, 2+ . . . , j+, . . . , n + ) , and ^(17, 2+ . . . , j+, . . . , n + ) . (2.23) 



III. REVIEW OF RECURSION RELATIONS 
A. On-Shell Recursion Relations for Trees 

The on-shell recursion relations rely on general properties of complex functions as well 
as factorization properties of scattering amplitudes. The proof [7] of the relations relies on 
a parameter-dependent shift of two of the external massless spinors, here labelled k and I, 
in an n-point process, 



where z is a complex number. The corresponding momenta (labeled by pi instead of ki in 
this section) are shifted as well, 




A/ + z\ k , 



(3.1) 





(3.2) 
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so that they remain massless, p\{z) = = pf(z), and overall momentum conservation is 
maintained. The shift also implies, 

h ^M z ) =t k - z (\n(k-\ + \k + )(i + \), 

f x ^p i {z)=p i + z{\r){k-\ + \k+){l+\). (3.3) 

Define a parameter-dependent continuation of an on-shell amplitude, 

A(z) = A(p u . . . ,p k (z),p k+1 , . . .,pi(z), . . . ,p n ), (3.4) 

evaluated at a particular set of complex momenta. When A is a tree amplitude or finite 
one-loop amplitude, A(z) is a rational function of z. The physical amplitude is given by 
A(0). 

Consider the contour integral, 

1 f dz 



where the contour is taken around the circle at infinity. If A (z) — > as z — > oo, as in the 
tree-level cases [6, 7, 8, 10, 31], then there is no "surface term"; that is, the integral (3.5) 
vanishes. Evaluating the integral as a sum of residues, we can then solve for A(0) to obtain, 

A(0) = - V Res ^ . (3.6) 

Z — Za Z 

poles a 

As explained in ref. [7], if A(z) only has simple poles, each residue is given by factorizing 
the shifted amplitude on the appropriate pole in momentum invariants, so that at tree level, 

^4(0) = J2 A l( z = z rs)j^A- h (z = z rs ) , (3.7) 

r,s,h r "' s 

where h = ±1 labels the helicity of the intermediate state. There is generically a double 
sum, labeled by r, s, over momentum poles, with legs k and I always appearing on opposite 
sides of the pole. The squared momentum associated with that pole, K^,.. s , is evaluated 
in the unshifted kinematics; whereas the on-shell amplitudes Al and Ar are evaluated in 
kinematics that have been shifted by eq. (3.1) with z = z rs , where 

K r-s 

(AT | $ r ... s \l~) ' 

To extend the approach to one loop [32], the sum (3.7) should also be taken over the two 
ways of assigning the loop to A L and A R . This formula assumes that there are no additional 
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poles present in the amplitude other than the standard poles for real momenta. At tree level 
it is possible to demonstrate the absence of additional poles, but at loop level it is not true. 

For the case of a fermionic pole, there is a sign subtlety similar to the situation with 
maximally- helicity- violating (MHV) vertices [4]. The correct fermionic propagator is, of 
course, ifc/k 2 . As is the case for MHV vertices, the $ is supplied by the amplitudes appearing 
on both sides of the pole. In these amplitudes, each momentum is directed outwards. Thus 
when we link two amplitudes across the pole we would obtain a numerator factor of the 
form |fci~)(&2~|, where k<i = —k\. This is not quite right since the same momentum argument 
should appear in both spinors, i.e., \k~l){k^\ ^ \kf ) ((— fci) + |. To correct this we flip the sign 
of the momentum in the spinor (A^~|. 

Because of the general structure of multiparticle factorization [39], only standard single 
poles in z arise from multiparticle channels, even at one loop. However, as was pointed 
out in ref. [32], double poles in z do arise at one loop due to collinear factorization. The 

splitting amplitudes with helicity configuration (+ + +) and ( ) (in an all-outgoing 

helicity convention) can lead to double poles in z, because their dependence on the spinor 
products takes the form [a b] / (ab) 2 for (+ + +), or its complex conjugate (a b) / [ab] 2 for 

( ) [38]. As discussed in ref. [32], this behavior alters the form of the recursion relation 

in an essential manner. In general, underneath the double pole sits an object of the form, 



which we call an "unreal pole" since there is no pole present when real momenta are used; it 
only appears, as a single pole, when we continue to complex momenta. As we shall discuss in 
section V, the finite quark amplitudes exhibit similar phenomena, except that in the quark 
case we shall also encounter unreal poles not directly associated with a double pole. 

B. One-loop Factorization Properties 

In order to build on-shell recursion relations, we need the factorization properties of one- 
loop amplitudes for complex momenta. It is useful to first review the factorization properties 
for real momenta, which we know from general arguments [39, 40, 46]. Although a good 
starting point, this is in general insufficient due to the appearance of unreal poles. 

As the real momenta of two external legs become collinear, any one-loop amplitude in 



[ab] 



(3.9) 
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massless gauge theory will factorize as, 

4 1} E (SplitL ^^,^; z) A^U ..,(« + 6)\ . . .) 

x=± 

+ Split?! (a^, b x »; z) A^. . . , (a + b)\ ...)), (3.10) 

where a and b are nearest neighbors in the cyclic ordering of legs. (When they are not nearest 
neighbors, there is no universal factorization behavior, but also no collinear singularity.) 
There are three distinct kinds of collinear limits to consider: two gluons becoming collinear; a 
gluon becoming collinear with a quark (or anti-quark); and a quark and anti-quark becoming 
collinear (in the case that they are adjacent). These limits simplify in the case at hand, 
where all external gluons have positive helicity, due to vanishing of relevant (n — l)-point 
tree amplitudes. The required splitting amplitudes are tabulated in ref. [38]. 

First consider the case when two cyclically adjacent gluons, a and b, become collinear. 
In this case behaves as, 

A L n {. . . , a + , b\ . . .) Splits (a+, 5+; z) A L n _,{. . . , (a + b)\ . . .) 

+ Split; ): > + , b+- z) A^U. . . , (a + b)~, . . .) , (3.11) 

where Split* 1 -*' 9 is the one-loop splitting amplitude with a gluon circulating in the loop. 
Similarly, A s n behaves as, 

A s n {. . . , a+, b + , . . .) -^U Split^V, b + - z) A s n _,{. . . , (a + &)+, . . .) 

+ Splitf*(a + , b + - z) A^(. . . , (a + 6)", . . .) , (3.12) 

where Split^ :s is the one-loop splitting amplitude with a scalar circulating in the loop. 
In both cases, the remaining two terms, with opposite intermediate-gluon helicity, vanish 
because Split +\a + , b + ; z) and A < ^1 1 (±, +,+,...,+) are zero. Now [38], 

Split£ ):ff (a + , b + ; z) = Split5 ):s (a + , b + ; z), (3.13) 

thanks to a supersymmetry Ward identity [33]. Taking the difference we see that in this 
limit the A^~ s class of amplitudes to be computed has a very simple structure, 

A L n ~ s (. . . , a + , b + , . . .) SplitL 0) (a + , b + ; z) A L n Z[{. . . , (a + b)\ . . .) . (3.14) 

In the limit that a gluon becomes collinear with the negative-helicity fermion, for either 
A 1 ^ or A s n we find, thanks to helicity conservation on a fermion line, that, 

A^(a 7 , b\ . . .) Splitj^- 6 + ; z) A L n ^((a + b)j , . . .) . (3.15) 
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Again taking the difference, we obtain, 

A L n ~ s {aj, b + , . . .) Split[° / ) )+ (a / -, b + ; z) A L n Z{{{a + b)j, . . .) . (3.16) 

A similar expression holds for the limit in which a gluon becomes collinear with the positive- 
helicity fermion. 

Finally, if the quark and anti-quark are adjacent, the amplitude factorizes onto the finite 
one-loop pure-glue amplitudes, 

A L n > s {aJ, bp . . .) Split^af, bp, z) A^a + b)\ . . .) 

+ SplitJ? (a~, bp, z) A^_ x {{a + bp , . . .) . (3.17) 

The gluon and scalar loop contributions to A^ 1 2 1 (±, +,+,...,+) are the same by supersym- 
metry [33]. Hence the difference A^~ s is finite in this limit, 

A^ s (aj, bp...) allb > non-singular . (3.18) 

In the case of multiparticle factorization, the vanishing of An\lq , 2+, 3 + , . . . , n + ) implies 
that we can only factorize on a gluon pole, 

A^(lf, T , r nr n + ) (3.19) 

4^(1//, 2 + , ..... lf /,,-. (- Kl ... m r)-^A^ m+1 (Kt m , (m + 1) + , . . . ,n + ) , 

A l...m 

with m > 3. Again A^~ s has a non-singular limit, 

A L n ~ s {lJ, 1 ./ f /,/• n + ) non-singular , (3.20) 

since once again the behavior of the L and s pieces are identical and cancel in the L — s 
difference. 

Because of the appearance of unreal poles (3.9) this is not the entire story. Unfortunately, 
as yet there are no general theorems to guide us on the factorization properties in this new 
class of poles. As shown in ref. [32] for the finite one-loop pure-glue amplitudes, factorization 
on the unreal poles is not solely into products of lower-point amplitudes; other factors arise. 
For the finite quark amplitudes we shall find, just as in ref. [32], a systematic set of correction 
factors. In section V we will comment on different types of factors that appear. 
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IV. REVIEW OF KNOWN FINITE QCD AMPLITUDES 



In this section, we collect previously-known results for tree and one-loop finite ampli- 
tudes, which feed into the recursive formulae for the finite quark amplitudes to be discussed 
in section V. 

We will need the tree-level MHV amplitudes [46, 49], 

4°> ( i+, . . . , mr , .... m2 -, ... , „♦) = i (12) ( ;^ )4 (nl) , (4.D 

and 

where we use the generic fermion label / again, and where the omitted labels refer to 
positive-helicity gluons. 

We will also need the one-loop pure-gluon amplitudes, either with all helicities positive, 
or with a single negative helicity. For the all-positive case with n > 4 legs, Chalmers and 
the authors [50, 51] wrote a conjecture [56] based on collinear limits, 

A® (1 + , 2+ . . . , n + ) = 1 - {12) (23> ...^_ 1)n) (nl) » ( 4 - 3 ) 



where 



and 



H n = - J2 Tr_[$JJJ i4 ], (A A) 



Tr. 



I<i\<t2<i3<ii<n 



= (M2) [kh] (kk) [kk] ■ (4.5) 



The conjecture was proven by Mahlon [37]. In eq. (4.3) we have extracted an overall factor 
of l/(47r) 2 compared with ref. [51], to be consistent with the normalization in eq. (2.2) for 
the quark amplitudes, and we have set a multiplicity-counting parameter iVp to 2, to match 
the number of color-stripped bosonic states defined to be circulating in the loop in the 
amplitudes A s n . More generally, the overall prefactor iVp is just the difference between the 
number of bosonic and fermionic states circulating in the loop. 

For n = 3 the proper vertex where legs 1 and 2 are the external legs is [32], 

4 , (1+ , 2+i3+) = iMMM. (4 . 6) 



K 2 
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It is useful to expose the kinematic pole, so we define the vertex, 



\/ 3 (1) (l + , 2+, 3 + ) ee K\ 2 Af (1+, 2+ 3 + ) = -i [1 2] [2 3] [3 1] 



(4.7) 



The single negative-helicity amplitudes may also be required as inputs into the quark 
recursion relation. The four-point amplitude was first calculated using string-based methods 
and is given by [52], 

(2 4) [2 4] 3 



4 1} (l-,2 + ,3 + ,4 + ) = 



(4.8) 



3 [12] (23) (34) [41] ' 

The five-point amplitude with a single negative-helicity leg was also first calculated using 
string-based methods and is given by [50, 53], 



4 1) (l-,2 + ,3 + ,4 + ,5 + ) 



[2 5f 



+ 



(14) 3 [4 5] (3 5) (13) 3 [3 2] (4 2) 



(4.9) 



3(34) 2 L [12] [51] ' (12)(23)(45) 2 (1 5) (5 4) (3 2) 2 . 
The n-point generalization was then obtained by Mahlon [37] via off-shell recursive methods. 
We recently used an on-shell recursion relation to obtain compact representations of the six- 
and seven-point amplitudes [32]. The six-point amplitude is, 



4 1} (l-,2 + ,3 + ,4 + ,5 + ,6 H 



(1"| (2 + 3) |6-)' : 



(i-| (3 + 4) \2-y 



(12) (23) (45) z s 123 (3"|(1 + 2)|6-) + (34) 2 (5 6) (61) s 234 (5"| (3 + 4) |2") 



+ 



[2 6f 



[2 3] [3 4] 



[4 5] [5 6] 



[3 5] 



[12] [61] s 3 45 \(4 5) (5"|(3 + 4)|2-) (3 4) (3"| (1 + 2) |6") (3 4) (4 5) 



(13) 3 [2 3] (2 4) 



(15) 3 (4 6) [5 6] 



(2 3) 2 (3 4) 2 (4 5) (5 6) (61) (12) (2 3) (3 4) (4 5) 2 (5 6) 2 
(14) 3 (35)(r|(2 + 3)|4~) 



(12) (2 3) (3 4) 2 (4 5) 2 (5 6) (61) 



(4.10) 



In section V, we give a compact expression for all n. 

The known results for the fermionic four-point amplitudes [54] are, 



^-(17,2/, 3+4+) = -- 
^(l / -,2 / + ,3 + ,4 + ) = 
and for the five-point fermionic amplitudes [35], 



(12) [2 4] 

2 (2 3) (3 4) ' 
j (lj) [2 4] 

3 (3 4) 2 ' 



(4.11) 



^-*(ly,2;,3 + ,4 + ,5 + ) 



% (12) [23] (31) + (14) [45] (51) 
2 (2 3) (3 4) (4 5) (51) 

17 



74*(l7,2+,3+,4+,5+) = -l( 



(13) [3 4] (4 l) 2 (14) (2 4) [4 5] (51) 



3 V (1 2) (3 4) 2 (4 5) (5 1) (1 2) (2 3) (3 4) (4 5)' 
[2 3] [2 5] 



+ 



7l£- s (17,2+,3+,4+,5- 



[12] (34) (45)7' 

1 <13> «12> [2 3] (31) + (14) [4 5] (51)) 

2 (12) (2 3) (3 4) (4 5) (51) 



74* ( l7,2+,3+4+,5+) = j (14)(15)[45] 2 . (4.12) 
5V /. » / > > ) 3 (12) (2 3) (4 5) 2 V ; 

V. QUARK RECURSION RELATIONS 

As was discussed in section II, for the finite helicity amplitudes — one quark pair, with 
all gluons having positive helicity — there are two independent primitive amplitudes that 
need to be computed, 

A L n - s {jf) ee 74^(17,2+,. .., J/ +,...,n+), 2<j< [(n + l)/2] , (5.1) 
A s n {jf) = 7^(17, 2+,.. ...,n + ), 2<j<n-2, (5.2) 

where we have adopted an abbreviated notation retaining only the label of the positive- 
helicity fermion. A computation of these primitive amplitudes then determines all the finite 
one-loop quark amplitudes in QCD. 

A. Structure of Five-Point Recursion for L — s Contribution 

Let us begin, as in ref . [32] , by examining the structure of five point amplitudes. Consider 
A^~ s {Zf), as defined in eq. (5.1), and choosing (k,l) = (1,5) as shift variables in eq. (3.1), 

Ai — > Ai , 

Ai — > Ai — zAs , 

A 5 -> A 5 + z\ 1 , (5.3) 
A5 A 5 , 

with all other spinors unchanged. Under this shift, using the known result for the amplitude 
we have 

a l- s(2+ . z) = (12) (13) [2 3] i (14) [4 5] 

5 {f,) 2 (2 3) (34) ((45) + ^ (41)) (15) 2 (2 3) (34) ((45) + ^ (4 1)) ' K ' J 
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FIG. 3: The real-pole diagram in the recursion relation for A\ s (l^, 2^, 3 + , 4 + , 5 + ). The vertices 
labeled by a T are trees, and the ones labeled by an L are loops. 

which vanishes as z — > oo, so no surface term is required. Looking at the z = pole in 
A^~ s (z)/z, we expect to find the following two terms, 

1 (12) (13) [23] i (14) [45] 

2 (23) (34) (45) (15) 2 (23) (34) (45)' 1 ' ' 
appearing in the recursion. 

On the other hand, from the structure of the collinear limits — eqs. (3.14), (3.16), and 
(3.18) — we expect to find only one term. A lone term, displayed in fig. 3, is indeed what 
emerges from the naive form of the recursion relation, 



D 1 = — ^(i / -,2;,3 + ,^ 4 + 5 )4 0, ((-K 45 )-,4 + ,5 + ) 
S45 

_ i 1 (12) [2 £45] [45] 3 

2 s 45 ( 2 3) (3 K45) [(-£45) 4] [5 (-K45)] 



% 1 (12)(1"| 


-^45 


2-) 


[4 5] 3 




2s 45 (2 3) (3"| 


-^45 


5->(l-| 


-^45 


4") 



(5.6) 

^ I -^45 I* / 
i 1 (12) (13) [3 2] [4 5] 2 
= 2^(2 3) (3 4) [4 5] (15) 
% (12) (13) [2 3] 
~ 2 (23) (34) (45) (15) ' 

exactly the first term in A^~ s . But the second term in eq. (5.5) is missing. 

What happened to it? Notice that the second term in eq. (5.5) is not singular in the fc 4 || k 5 
limit, so long as we are considering real momenta, because the numerator and denominator 
vanish at the same rate in the limit. If we consider complex momenta, however, the behavior 
of (4 5) is decoupled from that of [4 5], and in regions where (4 5) — > 0, there is a pole. 

In other words: while multiparticle and collinear factorization capture the full pole struc- 
ture for real momenta, they do not do so for complex momenta. In particular, a term 
containing the unreal pole (3.9) in the a \\ b limit, will not contribute to a collinear singu- 
larity, but will contribute to a pole for complex momenta (which decouple the behavior of 
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FIG. 4: The extra unreal pole contribution in the recursion relation for 2+, 3+, 4+, 5+). 

(a 6) from that of [a b}). That is, for our purposes we need to ask what the quasi-universal 
behavior in the collinear limit is of finite terms with non-trivial phase structure. For real 
momenta the ratio (3.9) is always nonsingular, because \(ab)\ = \[ab]\ = \J\2k a ■ (It 
does contain a phase dependence, though, which selects it out uniquely in the collinear limit: 

— * —* 

as the vectors k a and k b are rotated around each other by an angle the ratio (3.9) will 
change by the phase factor exp(— 2i<f>).) 

At tree level, a diagrammatic analysis that isolates the collinear singularities for real 
momenta extends readily (for massless amplitudes) to complex momenta. However, the 
same is not true at loop level. We will not offer such an analysis in this paper. Rather, we 
will exhibit a set of ansatze which use lower-point amplitudes and yield the contributions 
required for the corresponding poles. 

We proceed as in ref. [32]. Our experience there suggests that the missing term is asso- 
ciated with a factorization of the type shown in fig. 4. In order to define the 3-point loop 
vertex in this diagram, we could try to use the factor vj 1 ^ vertex (4.7) which appeared in the 
gluonic case [32]. However, eq. (5.5) contains a normalization factor of 1/2 rather than 1/3, 
and is not associated directly with a double pole. So we introduce the "L — s loop vertex", 

\/ 3 L - s (l + ,2 + ,3 + ) = -^[12][2 3][31] , (5.7) 

which is proportional to V^ l \ 

In the present case, there is no real-pole contribution associated with the diagram in 
fig. 4; it would correspond to an li A sn contribution which we have subtracted out in A L ~ S . 
(Accordingly, for complex momenta, the double pole contributions present in ref. [32] are 
absent.) However, unreal pole contributions can arise. There are a number of constraints 
that allow us to find the precise form of the required term. The unreal pole is power- 
suppressed by a factor of S45 compared to the diagram in fig. 4. By dimensional analysis we 
then need additional factors to obtain the correct overall dimensions. The additional factors 



20 



must be invariant under phase rotations of spinors associated with all external states, and 
the intermediate state K^. 

We can use universal multiplicative "soft factors," which describe the insertion of a soft 
gluon s between two hard partons a and b in a color-ordered amplitude, to construct such 
a multiplicative function. The soft factors depend only on the helicity of the soft gluon and 
are given by [46], 

S®(a,s+,b) = j^-jT, (5.8) 
{a s) {s b) 

SW(a,s-,b) = - J^tt. (5.9) 
[a s\ [s b\ 

They are invariant under phase rotations of spinors associated with a and b, but not s. 
However, the product 

s 45 <S (0) («, s + , b)S {0) (c, (s)-, d) (5. 10) 

is dimensionless, invariant under phase rotations of a, b, c, d and s, and suppressed as 
K\ h — > 0, for suitable choices of a, b, c, d and s. 

Since the leg s appears in two factors in eq. (5.10) carrying opposite helicity, it is natural 
to identify it with the on-shell intermediate momentum K 45 . Choosing c = 5 and d — 4 
produces the desired collinear behavior oc (4 5) for eq. (5.10). A little experimentation shows 
that here the legs 2 and 1 should be identified with a and b. Multiplying eq. (5.10) with 
these assignments by the "naive" diagram in fig. 4, 

4-4 0) (i/- 2+, 3 + , K^) V 3 L - S ((-K 45 ) + , 4+, 5 + ) , (5.11) 

gives us exactly the missing contribution, 

% (14) [4 5] 



(5.12) 



2 (2 3) (3 4) (4 5) 

More generally, in the n-point on-shell recursion relation for A^~ s (j^), when we choose 
the shifted legs in eq. (3.1) to be (k,l) = (l,n), we need to add an unreal-pole term of the 
form, 

1 (\- 9+ „•+ — o\+ ir- \m l ~s((_ £r. a+ („_i\+ ^ + 



^(l/, 2 + , . . . , j+, . . . , (n - 2) + , K( n _ 1)n )Vr s ((-K (n - 1)n ) + , (n - 1)\ h + ) 



V .C(0)(n £ 

>-l) 



x S^U, K+ iyn , 1) 5(°)(n, (-k {n -i)n)-, n - 1) . (5.13) 
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The complete recursion relation for this choice of shifted legs is then, 
A L n - s {3l) 

•//• (" 2) + , ^ { + n _ 1)n )4 0) ((-^(n-l)n)-, (n - 1) + , ^ 



^(n— l)n 



+ — — 4^1(17, (n- 2)+, xr_ 1) jy 3 L - s ((-x (n _ 1) „)+, (n - 1) + , h + ) 



s (n— l)n 

x S^a^.il^fM-^-iln)^-!), (5-14) 
where the hatted momenta in eq. (5.14) are defined by the shift 

Ai — > Ai — z\ n , 

A„ -> A„ + 2A1 , (5.15) 

with 

(r|l B -i,>-> «n-i)i>' 1 j 

in each term. This relation assumes that, as in the case n — 5, the shifted amplitude 
vanishes as 2; — > 00, so there is no surface term to be added to eq. (3.6). The recursion 
relation (5.14) can be solved, yielding the following simple expression, 

a l- s(1 - 2+ , + n+ x _ j M&n^Ain (517) 

A n {l f ,2 ,..., Jf ,...,n )-- (12) (23 >...(nl> ' ( j 

The reader may verify that this result has the correct collinear limits in all channels, and 
that it obeys the reflection symmetry (2.22). 



B. Structure of Recursion for s Contribution 

We turn next to the computation of the scalar contributions, given by the A s n function. 
In this case, there is a non-trivial collinear limit when k n _i \\ k n with real momenta, so that 
we expect to have both double-pole and single-pole contributions with complex momenta. 
The latter may have an interpretation as 'underlying' the double-pole term, or else as simply 
being unreal poles. The presence of both types of term is similar to the gluon amplitudes 
with one negative-helicity gluon considered in ref. [32]. We can again find an appropriate 
function using soft factors; the factors here are however different from those in the gluon 
case. As above, take the shifts to be (l,n). In the five-point case, for example, the factors 
are, 

^°)(3,X 4 + 5 ,2)<S(°)(5,(-X 45 )-,4), (5.18) 
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+ 

J f~ 
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(n-l)+ 



FIG. 5: The diagrams corresponding to the terms in the recursion relation in eq. (5.20). In diagram 
(a) I runs over {j + l,j + 2, . . . , n — 2}. Diagram (c) contains a double pole as well as an unreal 
pole underneath it. 

when the positive-helicity quark is leg 2, and absent if it is leg 3. This observation suggests 
a general form, 

5<°>(n - 2,k+_ 1)n ,j)S(°\n, (-£ (b _ 1)b )- n - 1) , (5.19) 

where j labels the positive-helicity quark. 

The full recursion relation for A^(jf~), again assuming the absence of a z — > oo surface 
term, is depicted in fig. 5, and reads, 



n-2 



= E ^A (0) (i/-, 2 + , ... ...,(/ - 1) + , K Ln )Ai:> l+2 ((-K L .. n )\ i\ ■ ■ ■ , ^ + ) 

+ —t—K-ifr, r , f (n 2)\ K+_ 1) JAf\(-K (n _ 1)n )-, (n - 1)+, h + ) 



|(!) 



S(n-l) 



+ -gJ—^i", 2' , /f (n 2) + , K^ 1)n )V^((-K (n „ 1)n )\ (n - 1)\ 



(n— l)n 



x (l + S(n _ 1)n S(°>(n - 2,£+_ 1)n , j)5^(n, (-^ (n _ 1)n )-, n - 1)) . 
The hatted legs undergo the shift in eq. (5.15), where in the sj... n channel we set, 

S/...n 

2 = 



(5.20) 



(5.21) 



(1"| #,...>-)" 

The tree-side soft factor here, — 2, -^^_ 1 ) n ,j), is in a sense the "complement" of 



that for s appearing in eq. (5.13); that is, 

5<°>(n - 2, j) + S%\ tf+ i) = 5<°>(n - 2, tf+ i) , (5.22) 



where the right-hand side is exactly the tree-side soft factor that appeared in our previous 
recursion relation [32] for the one-loop all-gluon amplitude with one negative helicity. 
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FIG. 6: The recursive diagrams for y4|(2^) using the shift in eq. (5.15) with n = 6. 
C. Solution to A s n Recursion Relation 

Unlike A^~ s , A s n does have multi-particle poles. Consider the six-point case. If we use the 
shift (5.15) the recursion relation (5.20) generates the diagrams shown in fig. 6 for the case 
where leg 2 is the positive helicity fermion leg. The other cases, where the positive-helicity 
fermion is leg 3 or 4, are similar. Evaluating these contributions, we find, 



t/ (14)(l-|(2 + 3)(3 + 4)|l+) (24) (15) (4 + 5)(5 + 6) |1+) 
3 V (1 2) (3 4) 2 (4 5) (5 6) (61) (1 2) (2 3) (3 4) (4 5) 2 (5 6) (6 1) 



(1 2) (2 3) (3 4) (4 5) (5 6) 2 (3 4) 2 (5 6) (6 1) (5" | (3 + 4) |2") 

[2 6] 2 (2+| (3 + 4) (4 + 5)(3 + 4) (4 + 5) |6~) 
[1 2] (3 4) (4 5) (5- 1 (3 + 4) |2") (3"| (4 + 5) |6") s 3 45 
| (2-|(4 + 5)|6-)(l-|(4 + 5)|6-) 2 \ 

(12) (23) (4 5) 2 (3-|(4 + 5)|6-) S4 56/ ' 
if (15) (16) (3 5) [5 6] (15) (Tj (4 + 5)(5 + 6) 

3 V (1 2) (2 3) (3 4) (4 5) (5 6) 2 + (1 2) (2 3) (4 5) 2 (5 6) (6 1) 



(25) (1"|56|1 + ) 



+ 



(1-|(3 + 4)|2-) 2 




(5.24) 



Aim 



i (1-|56|1+) 
3 (12) (2 3) (3 4) (5 6) 



(5.25) 



These expressions suggest the all-n forms, 



(n-l)n\l + ) 



(5.26) 



3(12) ■■■((n- 3) in -2)) ((n - l)n) 2 ' 
if (l-|(n-l)n|l + )((n-3)(n-l)) 
3 V (1 2) (2 3) • • • ((n - 2) (n - 1)) ((n - 1) n) 2 
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+ (l(n-l))(l-|^n-2,n-A-l,n|l + ) 

(1 2) (2 3) • • • ((n - 4) (n - 3)) ((n - 2) (n - l)) 2 ((n - 1) n) (n 1) 
+ { V,^T,) ln ~i, —2 ) ■ (5-27) 



(1 2) (2 3) • • • ((n - 4) (n - 3)) ((n - 2) (n - l)Y s (n _ 2) ... n . 

By studying the structure of eqs. (5.23)-(5.27), and also the output of the recursion 
relation for n = 7 and n = 8, we have arrived at the following compact formula, valid for all 

j and n, 

= l(12) S (23) S - 2 (nl) ' (5 ' 28) 



where 



51 " ^ (Uf~+i)> ' (5 ' 29) 

52 = <(/-!)/> 



, ; : „ , 1 t^y. n ti... v \(l ~ 1) + ) (1-| t^y-nti..., \V 

(P (p+1)) 



X 



and 



^ 2 ...(I-l)^i... p t2- { l-l)tl... v \(P + 1) + ) 

x(l~| Ki...p$(p + i)... n |l + > 2 (r| $i...p${p + l)... n |l + > 

?; (l-|^2-(t-l)[^,P)] 2 ^(H-l)-n|l + > (53Q) 
p-1 p 

^(Z,p)=5] ]T #J m . (5.31) 

j=Z m=i+l 

D. Verification of Solution 

We now verify analytically that these amplitudes satisfy the recursion relation (5.20). 
First consider the term shown in fig. 5(b), containing the three-point tree amplitude 
Af\(— K( n -i)n)~, {n — l) + )^ + )- Let Si and S 2 stand for the shifted versions of Si and 
S 2 for the appropriate (n— l)-point one-loop quark amplitude. The term in fig. 5(b) can be 
simplified to 

1 _ [(n-l)nf 

K l-i,n [nK n -i !n ] [K„-i, n (n - 1)] 

% s ± + s 2 

x- 



3 (ir n _ ljn l)(12) (2 3)-..((n-3)(n-2)) ((n - 2) K n ^ n ) 
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i St + S 2 ((n - 2) (n - 1)) ((to - 1) to) (to 1) [(to - 1) n] 3 

Z 3 (1 2) (2 3) • • • (to 1) ((to - 1) n) [(to - 1) n] (F\f^~ n \(n - 1)") <(n - 2)~| # n _ lifl \n~) 

i Si + s 2 

3(12) (23) •••(nl) ' 1 J 

We see that the correct spinor denominator factor is reproduced. 

The next question is how Si and S 2 are affected by the shift (5.15). Because Si is simpler, 
we discuss it first. Note that the (n — l)-point expression Si is a single sum over I containing 
(n — 2) — j terms, which is one fewer than the number of terms in the n-point expression Si 
we are trying to produce. All terms but the last in Si have very simple behavior under the 
shift (5.15) of Ai and A n . They depend on Ai through (1~| and |1 + ), but do not depend on 
Ai. Their dependence on A n is solely via the factor 

• • • (#( I+ i)... (n _2) + tn-l,n)|l + ) = • • • ^(m)...n|l + ) , (5-33) 

because the shift in ^t n _ l n is proportional to A x . Thus each such term directly yields the 
corresponding term in the n-point sum Si. 

The last term in Si, with / = n — 2, is the exception, because it has additional dependence 
on ^t n _ l n . It can be written as 

(j(n-2)) (lK n .l, n ) [ATn-l,nn](l-|(4-2 + ^n-l,n)^n-l,n|l + ) 

((n - 2) K n - ljn ) [K n _ hn n\ 
( j(n _ 2 ))(l(n-l))(l-|4_ 2 ^ ra . liTO |l+) 
((n-2)(n-l)> 

which does not quite match the next-to-last term (I — n — 2) in Si, 

(j (U - 2)) (1 (TO - 1)) # B - 2 ,n-l#n-l,n + ) 

((„-2)(n-l)> 

Additional contributions to the next-to-last term, and the whole of the last term (I = 
to — 1) in Si, come from the term containing V^ 1 \(—K^ n _i^ )n ) + , (to — 1) + , to + ) in the recursion 
relation (5.20), depicted in fig. 5(c). This term can be simplified to, 

2 

j 1 (1 ifn-l.n) (j Kn-lfid [-^n-l.n (w ~ 1)] [fo ~ 1) wj [n K n -i, n ] 

3 s* n _ 1)n (1 2) (2 3) • • • ((to - 3) (to - 2)) ((to - 2) K n _ hn ) 

((to-2)j)[to(to-1)] 



(5.34) 



(5.35) 



X 1 + S („-!)„ 



((to - 2) #„_!,„) (K n _i, n j) [nif„_i >n ] [if„-i,„ (n - 1)] 
1 



3(12) (2 3) •••((n- 3) (to - 2)) 
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((n-2)"| 


$n-l,n 


n~) 0'~l 




n->(l-| 


^n-l,nl(n-l)-> 



Z 1 (7(7l-l))(l7l)(l-|^ w _ 1|n 4|l + ) 

3(12) (2 3) ••■(nl) ((n-l)n) 
(j(n-2))((n-l)n)(l(n-l)) 



x 1 + 



(;(„- 1)> ((„- 2) (n-l)>(ln>. 
1 /0-(n-l))(ln)(l-|^ n _ 1 X|l + ) 



3(12) (2 3) •••(nl) V ((n-l)n) 

(j(n-2))(l(n-l))(l-|^ n _ 1 ^ n _ lin |l + ) 



(5.36) 



(( n _ 2 ) ( n -l)> 

The term containing (j (n — 2)) in eq. (5.36) combines with eq. (5.34) to generate the next- 
to-last term (5.35) in Si, while the term containing (j (n — 1)) is just the last term in Si. 
So we have demonstrated that all of the terms in Si are produced correctly by the recursion 
relation. 

Now we turn to £2- We divide the S2 terms into those with p < n — 1 and those with 
p = n — 1. The bulk of the terms, with p < n — 1, come from the S 2 contribution we 
left unexamined in eq. (5.32). To show this, we observe that (as was the case for the bulk 
of the Si terms) Ai never appears in S 2 - Also, k n only appears in eq. (5.30) for S 2 via 
(l _ |^t( p+1 )... n . . . or . . . ^t( p+1 )... n |l + ). Thus we may again apply eq. (5.33) (with / replaced 
by p) to S2, in order to see that every term with p < n — 1 in S 2 is generated directly from 
the corresponding term in S 2 . 

The S 2 terms with p = n — 1 come from the terms containing the one-loop pure-glue all- 
plus amplitude A n _ l+2 in eq. (5.20), depicted in fig. 5(a). The I th term in S 2 comes from the 
/ th term in the recursion relation. The main technical detail is to rewrite the numerator H n 
of the all-plus amplitude, as given in eq. (4.4), so that it looks like the operator [IF {I, n — l)] 2 
appearing in the terms with p = n — 1 in S 2 . For this purpose, we first use the identity (A12), 
with p — n, to rewrite the numerator factor H n -i +2 (l,n) (for the all-plus amplitude with 
gluon momenta ki, ki +i , . . . , k n , P, where P = —Ki... n is on-shell) as, 

H n . l+2 (l,n) = -lTr_[^(/,n)] 2 . (5.37) 



Then we use the identity 

2 L ~ " /J n K,.. n 
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^_ [Hl , n)f= ^m^m^, (5.38) 



where T is defined in eq. (A6), which follows similarly from eqs. (All) and (A5) for p = n. 
Finally, eq. (A20) shows that we can shift n — > n — 1 in the argument of T in eq. (5.38). 
Combining this equation with eq. (5.37), reversing the order of the spinor strings, using 
momentum conservation, and multiplying numerator and denominator by common factors, 
we arrive at the identity, 

H n _ l+2 (l, n) = —— — — . (5.39) 

(1 \$2-{l-l)K\^ + ) 

A few other identities are also useful, 



in 



1) Sl ... n = # 2 ...(i-l)#i...(n-l) K) , (5.40) 



(n-\^... n \n-) = (n-|#,... w |n-)- n _.^" n . rr (l" | ti-n \n~) 



\n 



-Si-..(n-l) , (5-41) 



<(n-l)n) = ((n-l)n)- ^ 2 "» ((n - 1) 1) 

(1 \Wi... n \n ) 

(l-|^ 2 ...(i-i)^i...(„-i)|(n-l) + > 



(1 |#j... n |n ) 

Using these identities, the I th term in the recursion relation can be written as, 

(I KL.. n ) 2 (j K t ... n ) 

3 *,...„ (12) (2 3) •••((/ -2) (J-l)) ((Z-l) £,...„) 

H n _ l+2 (l,n) 



(5.42) 



x 



(#,...„ Z) (/(/ + !))■ ..((n-l)n) (n 
z 1 ((I- 1) Q ((n - 1) n) (n 1) \n-f |n~) 

3 (1 2) (2 3) • • • (n 1) (n + | ^,.. n |(Z - 1) + ) (n + | ^,.. n |Z + ) 

(nl) (l-|^,... n |n-) 



x 



X 



(1 I #2-(I-l)#/-(n-l) l™ + ) (1 I #2...(!-l)#J...(n-l) I (" ~ + ) 

1 (l-|^ 2 ...(,_ 1) [^(/,n-l)] 2 ^|l+) 

«/-(n-l) (l~|-^2...(i-l)4|l + ) 

i ((*-!)*) 

3 (1 2) (23) • • • (n 1) (1"| 4^...(n-i) |(Z - 1) + ) (1"| *„#,...(„-!) IO 

{{n — 1) n) 



x 



(l-|^2...(/-l)^...(n-l) \{n~ 1) + ) (l-|^2...(i-l)^ I ...(n-l) l™ + > 
X (1-| tl.. <n -4n |1 + ) 2 (r|^...(n-l)^ + ) 

x (11^ 2 ... a - 1 )[^(/,n-l)] 2 4|l + ) (5 
Sj-(n-l) 

The final form is just the Z th term with p — n — 1 in S 2 in the solution (5.28). This completes 
the demonstration that eq. (5.28) obeys the recursion relation (5.20). 
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E. Factorization Properties of Solution 

In appendix A we show that formula (5.28) has all the correct multi-particle poles, fac- 
torizing properly onto products of the quark-containing MHV tree amplitudes (4.2) and the 
one-loop all-plus-helicity pure-gluon amplitudes (4.3). 

The collinear singularities of eq. (5.28) are also quite manifest, although we shall not 
verify them all in detail here. Most of the factors that diverge in the collinear limits are 
contained in the 1/((1 2) • • • (n 1)) prefactor. The generic collinear limit factorizes onto 
another quark amplitude in the same sequence (5.28), but with one fewer external gluon. 
If the two particles becoming collinear are color- adjacent gluons labeled I and I + 1, with 
I > j, then the amplitude can also factorize onto the product of the helicity-flip loop splitting 
amplitude, Split (1+ 1)+; z) cx [I (I + 1)] /(/ (/ + l)) 2 , and a quark-containing MHV tree 
amplitude (4.2). The corresponding terms in A s n (jj~) come partly from the / th term in S\, 
and partly from the term in S% with p — l + 1. The latter term has an apparent singularity 
from a manifest factor of However, this factor is cancelled by the numerator factor 

containing + l)] 2 = $i$i+iki$i+i = x + 1). On the other hand, the 

denominator factor (1~| W2-(i-i)Wi- P \p + ) produces a spinor product (/ (/ + 1)), like the one 
manifest in the I th term in Si. 

When the positive-helicity fermion is in the position j = 2, the limit 1 || 2 (where the anti- 
quark and quark momenta become collinear) factorizes the finite quark amplitudes (5.28) 
onto the finite pure-glue amplitudes, according to eq. (3.17). A term containing 1/ [12], 
corresponding to the factor Split^(l^, 2^~; z) in eq. (3.17), will factorize onto the all-plus 
amplitudes (4.3). The relevant term in A s n (2^) is the term in S2 with I = 3 and p — n — 1. 
The denominator factor of [1 2] comes from the string (1~ | W2~-(i-i)Wi--- P I ip + 1) + ) appearing 
in this term. This limit can be checked by analysis very analogous to the discussion of the 
multi-particle factorization limit Kf,,, n — > at the end of appendix A. 

The terms containing a factor of 1/ (12), corresponding to the factor Split+^l^, 2^; z) 
in eq. (3.17), will factorize onto all-gluon amplitudes with one negative-helicity gluon, 
A { n\l-,2 + ,...,n + ). We can use this factorization to extract a compact form for these 
pure-glue amplitudes. In this case, virtually every term in Si and S2 contributes, except for 
those terms in S 2 with I = 3. (In each of the / = 3 terms, the factor of 1/ (12) from the 
prefactor is cancelled by factors in the S 2 term itself.) The compact form for the n-gluon 



29 



amplitude that we derive from this limit is very similar to the quark amplitude itself (with 
3 ~* 1): 

4"(l-.2+,3V.., n +)4 (ia g + \ n . (5-44) 



where 



rp ^ (lQ(l(/+l))(l-|^, m ^ +1 )...Jl + ) 

1=2 

n— 2 n—1 



(1(1 + 1)} 

((1-1)1) 

t, ... It, \( 

1=3 p=l+l 



X 



ft(p+l)-nftl-p \(l ~ !) + ) (1 I ft(p+l)-nftl-p \ l ~ 
(p (p+1)) 



\-\t2- { l-X)tl... v \P + ) (l-\$2...(l-l)tl... p \(P + + ) 

x(l-|^...^ (p+1) ... n |l + > 3 

Jl-|^ 2 ...(^)[-F(^)] 2 ^( P+ l)...n|l + ) ^ {5M) 
Sl-p 

This result agrees with eqs. (4.8), (4.9) and (4.10), although it is written in a slightly 
different form. Furthermore, we have checked numerically that it agrees with the all-n 
result of Mahlon [37] up through n — 18. 



VI. FURTHER CROSS CHECKS 



There are a number of non-trivial checks we have performed on our results. As we already 
discussed, they have the correct factorization properties in all collinear and multiparticle 
factorization channels (with real momenta). Another powerful check, which we describe 
now, is that the amplitudes can be used to obtain formulae for certain known results for 
QED and mixed QED/QCD amplitudes, and they agree with those earlier results. 

We also comment on the consistency of using various shift variables. In particular, there 
are choices of shift without the subtlety of unreal poles, which lead to alternative recursion 
relations. These relations are also satisfied by the all-n expressions (5.17) and (5.28). 



A. Checks Based on QED Amplitudes 

Mahlon [36, 37] has computed the one-loop amplitudes for two separate processes, which 
can be related to the ones presented here by converting gluons into photons using appro- 
priate permutation sums. These results therefore provide a stringent cross check which is 
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independent of the subtleties associated with unreal poles. A more general discussion of how 
to convert primitive QCD amplitudes with an external qq pair into QED amplitudes may 
be found in appendix D of ref. [35]. 

In ref. [36], the pure QED amplitudes for e + e~ — > 7+7+ ... 7+ were computed. In ref. [37], 
the mixed QED/QCD amplitudes for e + e~ — > g + g + ■■■g + , and for r ) ± g + g + . . . g + , via a 
massless quark loop, were presented. Both sets of computations were performed using a 
recursively-constructed tree- level current for two off- shell massless fermions and an arbitrary 
number of gluons (or photons) [37, 55]. To obtain the one- loop amplitudes, this current is 
"sewn up" into a loop by joining the two off-shell ends with the vertex for emission of a real 
photon, or of a virtual photon coupled to an electron-positron pair. 

Consider first the case of e + e~ — > 7+7+ ... 7+. Here there are two types of contributions, 
termed Ai and A2 in ref. [36]. The A± contributions include a closed fermion loop, while 
the A2 contributions do not. Both contributions allow for photons to be emitted off the 
external fermion line, as well as from the closed fermion loop (in the case of A\). The A\ 
piece is related to our A s amplitude, while the A2 piece is related to A L ~ S . 

The A2 piece is a bit simpler and can be checked analytically, so we begin with it. 
Consider the primitive QCD amplitude A^(n^) = A^(lj, 2 + , . . . , (n — 1) + , rif). According 
to eq. (2.15), there is no fermion or scalar (or gluon loop) contribution in this case (simply 
because there are no external gluons on the same side of the external fermion line as the 
putative closed loop, and tadpole diagrams vanish here). Hence A s n (n^) = 0, and A^(n^) = 
^n _s ( n /) i s gi ven by eq. (5.17) for j = n. To convert this primitive amplitude into the 
QED amplitude A2 with no closed fermion loop, we merely need to sum eq. (5.17) for j — n 
over all (n — 2)! permutations of the n — 2 gluons. The sum over permutations cancels the 
diagrams with gluon self-interactions, but retains abelian emission off the "left" side of the 
fermion line. Thus we have, 

i (V2er V (In) gg g^^ffl 
" 2 (4tt)2 (1 a(2)> (a(2) a(3)) • ■ ■ (a(n - 1) n) (n 1) ' 1 " j 

UfcO„_2 

where the sum over a runs over all permutations of legs 2,3, ... ,n — 1, the subscript e 
signifies an electron, and the subscript 7 signifies a photon. We have replaced the QCD 
coupling g by the QED coupling \/2e, where the extra is due to our normalization of the 
color matrices ( Tr (T a T b ) = 5 ab ). This result matches A2 as given in eq. (77) of ref. [36], up 
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to an overall factor of (— l) n+1 . A factor of (— l) n in this difference probably comes from an 
opposite sign convention for the polarization vector of a positive-helicity massless photon 
(or gluon). The remaining sign may arise from an external-fermion sign convention. (We 
use a convention compatible with supersymmetry, as described in ref. [35].) 

Now consider Mahlon's Ai piece having closed fermion loops. The conversion of our 
primitive amplitudes with a closed fermion loop to a QED amplitude is again given by a 
permutation sum over all gluon legs to convert them to photons, 

A QED, fermion loop^^ 2 + . . . , (n - 1)+ 71+) 

= J -^f- E K(l F ,n+,a(2)\...,a(n-l)+), (6.2) 

where the permutations er are the same as in eq. (6.1). The overall sign accounts for the 
sign difference between a scalar and fermion in the loop. The factor of (—1)™ takes into 
account a minus sign between gluons emitted on the left side of the fermion line, and those 
emitted on the right side. Eq. (68) of ref. [36] gives A\. If we multiply that result by the 
same (— l) n+1 normalization factor as for A2, and convert it to our notation, we obtain, 

^QED, fermion loop^- 2 + . . . , (n - 1)+ 71+) 

= 3 (4tt)' ^ (ct(2) ct(3)) (ct(3) ct(4)) (ct(4) ct(2)> (1 a(5)> (cr(5) cr(6)) • • • (a(n - 1) n) 

X (K(2) + k a (3) + K{±)) 2 ' 
We have confirmed numerically that this result matches eq. (6.2) through n — 9. 

We have also recovered the same QED amplitude (6.3) using primitive amplitudes with 

the positive-helicity fermion at other locations in the ordering, by adjusting the combinatoric 

factors appropriately. For the class of A s n primitive amplitudes where the positive-helicity 

fermion is in the j th position, as in eq. (5.28), we first relabel the amplitude so that the 

gluons still run from 2 to n — 1, by letting j—*n,j + l—*j, etc. Then we sum over 

permutations according to, 

^QED, fermion loop (l - 2 + . . . , (n - l)+ 71+) 

= ~ ( ~y E ^(l7^(2) + ,...^(j-l) + ,n + ,cT(j) + ,...,cT(n-l) + ).(6.4) 

The sign factor is (— 1) J because now only n — 2 —j gluons are emitted from the "right" side 
of the fermion line. 
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The combinatoric factor ("Zg) can ^ e ex pl &me d simply if we assume that all QED con- 
tributions vanish whenever more than four photons attach to the massless electron loop - 
one photon virtual, the rest real and of positive helicity. This assumption is quite reasonable 
because the multi-photon amplitudes ^4 QED (7 ± 7 + • • -7+) vanish for all n > 4, for either sign 
of the first photon's helicity [36]. The contributions with fewer than four photons vanish by 
Furry's theorem, and the vanishing of massless external leg corrections. Then of the n — 2 
external gluons in the primitive amplitude, precisely three must attach to the fermion loop, 
and n — 5 to the external fermion line. The latter n — 5 gluons are divided into j — 2 on 
the left side and (n — 5) — (j — 2) = n — j — 3 on the right. When we sum over all gluon 
permutations, we overcount QED diagrams that differ only by a re-ordering of photons on 
the left with respect to those on the right, which preserves the order within the left set, and 
within the right set. This overcount is (n — 5)1/ (j — 2)1/ (n — j — 3)! = ("^ )• Conversely, we 
can take the fact that eq. (6.4) works numerically (for all j and n < 9) as evidence in favor 
of the vanishing of the one-off-shell, (n — l)-positive- helicity photon amplitudes for n > 4. 

We also may use our primitive amplitudes to construct mixed QED/QCD amplitudes for 
an electron-positron pair plus (n — 2) positive-helicity gluons, and compare with Mahlon's 
earlier computation [37]. In this case, the permutation sums are a bit more involved. The 
mixed amplitudes are obtained using two separate sums, which cancel contributions where 
a gluon is attached to the electron line, and which allow for all possible orderings of the 
virtual photon attaching the e + e~ pair to the loop, with respect to the color-ordered gluons. 
We again let the negative- helicity fermion be leg 1, but label the positive-helicity one by 2 
(instead of n). 

In the five-point case, for example, the appropriate sum for the coefficient of the color 
trace Tr(T a3 T a4 T as ) in ^5(17, 2+ 3 + , 4 + , 5+) is, 

A— d (l e -,2 e +,3+,4+,5+) = -^^^(^(17,2;, ,7(3)+ a(4)V(5)+) 

+ ^(17, a(3) + , 2/, a(4) + , a(5) + )) , (6.5) 

where the sum runs over the cyclic permutations, Z 3 , of the gluons legs {3,4,5}. The 
unwanted diagrams appear in pairs in the permutation sum, but with opposite signs due to 
the antisymmetry of color-ordered vertices. 

More generally, the coefficient of the color structure Tr(T a3 T a4 . . .T a ") in the n-point 
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amplitude A n (l e , 2+, 3 + , . . . , n + ) is, 
<r d (l e " 2+ 3+,...,n + ) 

/ /q \2 n— 2 n— 2 

W - rn/ " E E -(3) + , • • • , -(J) + , 2+ <r(j + 1)+, . . . , a{nY) . (6.6) 



(4tt) 2 

The permutation sum is over cyclic permutations, Z n _ 2) of the n — 2 gluon legs, labeled by 
{3,4, . . . , n}. The sum over j is over the primitive amplitudes A s n {j^) with the positive- 
helicity fermion in the j th position. (The primitive amplitudes with fewer than two trailing 
gluons, j = n — 1 or n, vanish trivially and have been dropped from the sum.) 

Although Mahlon's corresponding formula, eq. (53) of ref. [37], contains some errors, it 
is not difficult to use his eqs. (27) or (33) to rederive an expression for his form of the 
amplitude. (The limit on the final sum in eq. (33) should start at £ = 2, not i — 1.) With 
our normalization conventions and labelling, the result reads, 

n;1 ^ e ' e ' 3 ( 47r )2 (34)...((n-l)n)(n3) Sl2 

x E [M3)](-(°&)-\*ntm\ 1+ ) ( 6 - 7 ) 

n-1 

+ 2^(a(3)-|t((( + i)...n) ^(0l 1+ >)- 

The permutation sum again runs over cyclic permutations of the n — 2 gluon legs, 
{3,4, ...,n}. We have confirmed numerically through n — 16 that our expression (6.6) 
agrees with eq. (6.7). This check is particularly useful because all primitive amplitudes 
A^(jf~) enter into the permutation sum. 



B. Consistency of Various Shifts 

We can also make different choices for the shifted momenta (3.1) used to derive the 
recursion relations. For many choices, the unreal poles found in the previous section lead to 
extra terms with different soft-factor coefficients, which we again determined empirically. In 
other cases, the unreal poles are absent entirely. These latter choices give us an independent 
check on the factors described above, because there are no correction factors to be determined 
at all. Only the "naive" terms, with the expected single-power denominators, enter. It turns 
out that we can make such "clean" choices for all amplitudes A s n with j > 2. 
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FIG. 7: The two recursive diagrams for A s n with j > 2, using the shift in eq. (6.8). Diagram 
(b) vanishes. Unreal poles do not appear in diagram (a). 

An especially nice choice with this property is the (k, I) = (3, 2) shift, 

A3 — > A3 — 2A2 , 

A 2 -> \ 2 + z\ 3 . (6.8) 

With this shift only the two diagrams displayed in fig. 7 are generated; other potential dia- 
grams do not contribute because the tree amplitude appearing in the factorization vanishes 
or because they would contribute to R type primitive amplitudes instead of the L type 
under consideration. Since the diagrams factorize onto tree-level (not one-loop) three-point 
vertices, they do not contain any unreal poles. Very neatly, diagram (b) vanishes because 
the three-point tree Af\(—K 3A )~ , 3 + , 4 + ) vanishes with the shift (6.8), leaving only diagram 
(a). Furthermore, the evaluation of diagram (a) is rather simple, being essentially the same 
as for MHV tree amplitudes. For A s n) inserting the value of the (n — l)-point amplitude 
from eq. (5.28) into this diagram immediately yields the corresponding n-point amplitude, 
providing a simple verification of eq. (5.28) for j > 2, with the j = 2 boundary case as input. 

For j = 2 the shift (6.8) generates a surface term, since the shifted amplitude does not 
vanish as z — > 00. However, the surface term can be taken into account easily for checking 
a result, if not deriving it, as was demonstrated for the all-plus n-gluon amplitudes [32]. 
There is no unreal-pole contribution because the vertex V^^(3 + , 4 + , (— K^) + ) vanishes in 
the shifted kinematics. We have checked that this recursion relation is satisfied by the 
solution (5.28) for j = 2. 

Similarly, for the case of A^~ s (j^), with the same shift, (k,l) = (3,2), we also can 
construct a recursion relation with a surface term, which is obeyed by eq. (5.17) for j > 2, 
again with input from the j — 2 case. (For j = 2, this shift produces an A(z) which diverges 
as z — > 00, so no useful recursion relation is obtained.) 
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Numerous other choices of shifts are possible. For many, if not all, choices of shift 
variables, the contribution of an unreal pole will have the form (in the s a b channel), 

± — A R V 3 (a + , b + , (-K ab ) + ) S<°> (b, (-K ab y, a) (V a S<°> (u u K+, Vi j) , (6.9) 

i 

where the coefficients q are either 1 or 2, Ar is a tree amplitude and V3 is the three-positive 
helicity loop vertex. (We can make the relative signs positive by interchanging ui <-> Vi if 
necessary.) 

In summary, in all shifts that we have checked, one can obtain consistent results by 
appropriately adjusting the unreal pole contributions. For j > 2 we avoid the unreal poles 
with a suitable choice of shifts. More stringent tests along these lines would require a first 
principles derivation of the factors appearing in unreal poles. We defer this to future study. 

VII. CONCLUSIONS 

In this paper, we have presented compact formulae for QCD amplitudes with one quark 
pair, and n — 2 gluons of identical helicity. As explained in section II, these amplitudes 
are built out of primitive amplitudes A^ s and A s n . Our principal results are eqs. (5.17) 
and (5.28), for the all-n forms of these primitive amplitudes with positive-helicity gluons. 
(The corresponding primitive amplitudes with negative-helicity gluons can be obtained by 
parity, implemented by spinor conjugation.) We also provide a new compact representa- 
tion (5.44) of the previously-obtained [37] n-gluon amplitudes with a single negative helicity 
and the rest positive. 

The corresponding tree-level quark-gluon amplitudes vanish, and hence these amplitudes 
are both infrared- and ultraviolet-finite. These were the last unknown finite loop ampli- 
tudes: formulae for the finite n-gluon amplitudes [37, 51] have been known for a while, 
and amplitudes with additional quark pairs or higher loops are necessarily divergent in four 
dimensions. 

Because the corresponding tree-level amplitudes vanish, the finite one-loop amplitudes we 
computed do not enter next-to-leading order cross sections for jet production at hadron col- 
liders. However, they do contribute at next-to-next-to-leading order. The finite amplitudes 
also appear in factorization limits of the remaining, divergent one-loop QCD amplitudes. 
Accordingly, their structure will likely play a role in understanding the latter amplitudes. 
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We constructed these amplitudes via on-shell recursion relations. The construction of 
such relations relies on knowledge of the factorization properties of amplitudes in complex 
momenta. At tree level, these properties are determined by the factorization properties 
in real momenta, which are known to be universal. At loop level, this is no longer true. 
Generic loop-level relations differ from tree-level ones in having "unreal poles" - - poles 
that are present for complex momenta, but absent for real momenta. In lieu of analogs 
of the standard real-momentum factorization arguments [38, 39, 40] for this class of poles, 
we took a heuristic approach. We empirically determined the structure of terms in the 
recursion relations associated with the unreal poles with the aid of known four- and five- 
point amplitudes. Then we applied the same structure to the case of additional external 
legs. For most (but not all) of the primitive amplitudes, we were able to find choices of 
complex shifted momenta which avoid the unreal poles. The agreement of the results with 
these alternate shifts provides a strong consistency check on the approach we took. In 
addition, we verified that our results satisfy the required collinear and multiparticle-pole 
factorization forms (in real momenta). A third and independent stringent check comes from 
a comparison of certain QED and mixed QCD/QED amplitudes computed by an entirely 
different method. 

Unreal poles are an essential feature of the analytic structure of loop amplitudes which 
deserves further study. A first-principles understanding of the extent of their universality 
and the structure of factorization would be very important. Such an understanding would 
strengthen the use of loop-level recursion relations as a complement to the unitarity-based 
method for performing loop calculations. 
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APPENDIX A: MULTI- PARTICLE FACTORIZATION OF A s n (j+) 

In this appendix we verify that the amplitudes (5.28) have the correct multi-particle 
factorization properties. As mentioned in section III B, these helicity amplitudes only have 
multi-particle poles with gluonic intermediate states, not fermionic ones. Given that the 
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sum S 2 in eq. (5.30) contains manifest factors of it is convenient to first consider the 

limit K?... p — > where j < I and I + 1 < p < n. All the non-trivial multi-particle poles are 
covered by this case, except for those with p = n, which we shall discuss subsequently. 
As Kl„ p — > 0, we expect to find that 

Kijf ) ^ 1 •//• !) + > (~ P r, (p + • • • ,n+) 

x-^A« +2 (P + ,Z + ,...,p + ), (Al) 

where P = —Ki... p , and A p ^ l+2 is the all-plus pure-glue amplitude (4.3). Let H p _i +2 (l,p) 
denote the numerator factor H n defined in eq. (4.4), after the external momenta are relabeled 
to correspond to A p 1 }_ l+2 (P + , l + , . . . ,p + ). Using also eq. (4.2) for the quark-containing tree 
amplitudes, we expect the behavior 

K.^O i (iPf(jP) 



-> — - 



3<12>--.<(Z-2) (I - 1)> ((I - 1) P) (P(p + l)>...(nl> 

X J H p . l+2 {l,p) 

si... P (PI) (1(1 + !))■■■ ((p-l)p)(pP) ' 1 ] 

Now examine the term labeled by / and p in eq. (5.30) for S 2 , in the multi-particle factor- 
ization limit. Using K 2 ...(i-\) = —ki + P — K( p+ iy.. n , we have 



S2(l, P ) 



((1-1)1) 



(1"| t {p+1) ... n \P~) (P (I - 1)) (11 # {p+1 y.. n \P~) (PI) 
x (P(P+1)) 

$ {p+ ,y.. n \P~) (Pp) (1-| $ {p+ ry..n \P~) (P (P + 1)) 

x(l P)\l-\$ [p+1) ... n \P'f (jP) (l-\$ [p+1) ... n \P~) 

JiP)(i-\$ {P+1) ... n [Hi,p)\ 2 \P~) (A3) 

s l---p 

where we also used the fact that 

(X + \[^(l,p)] 2 \X-)=0 (A4) 

for any spinor, or spinor string, X. Equation (A4) follows from a more general result, 

<X+| [Hl,p)f \Y~) = - (Y+\ [Hl,p)f \X~) = - (Y + \ [Hl,p)f \X~) , (A5) 

where the lengths of strings X, Y are the same mod 2 (otherwise there is an additional sign 
in reversing them). Here we have defined the reversal of jF(/,p), 

p(i, P ) = P ^2J2 (A6) 

3=1 i=j+l 
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so that 



V 2 
i=l 



(A7) 



in the factorization limit. Comparing the limiting behavior (A3) of S 2 with the expecta- 
tion (A2), and cancelling various factors of W(j>+i)- n we see that the limit will be 
correct if we can show that 



ft( P +l)-nlHl,P)} 2 \P~) = ~ t^+iy-n \P-)H p _ l+2 (l, P ) . 

The Schouten identity, 

[a b] [cd] = [ac] [bd] + [ad] [cb] , 
together with eq. (A5), imply that 

(r\^ iP+1) ... n [T(i,p)] 2 \p-) = (r\^ iP+1) ... n \p-)Tr4T(i,p)] 2 

-(r\# (jH . 1) ... n [r(i, P )] 2 \p-) , 



(A8) 



(A9) 



(A10) 



or 



[l-\l£ (p+1> .. n [F(l, P )} 2 \P-) = - {l-\$ {p+1) ... n \P-) Tr_[^(/,p)] 2 . (All) 



So we just need to show that 



Tr4F(l,p)] 2 = -H p „ l+2 (l,p). 



(A12) 



We first use momentum conservation to remove the terms in H p -i +2 (l,p) (see eq. (4.4)) 
which contain the massless leg P: 

p 



Hp-i +2 (l,p) E Tr 

l<ii<i2<ia<i4,<P 

i 3 -l 



i- it i- ll- — \ \ t> ]L i. it. i. 

Nl fW2 N3 N4 / , / j 11 - ^i 1 r^i 2 ^i3^i4 

I<il<l2<l3<p 14=1 

= - E E^-fcAWj- (A13) 

I<il<t2<i3<p 14=1 

On the other hand, the left-hand side of eq. (A12) can be rewritten using the Schouten 
identity (A9) as, 

l - Tr_ [F(l, V)f = \ E E M fa <s] (<s <4> [U k] 

l<ii<i2<p l<i:i<i4<P 

= \ E E ( _ ^ ^ ^] ^ 4 ^ ^ + fa ^ fa ^ fa ^ fa ^0 

(<il<»2<P I<t3<l4<p 

= ~\ E E (M2) [M3] (M4) [Mi] ■ (A14) 



l<il<t2<P (<«4<*3<P 
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The last term in the second line vanishes using Kf = 0. On the right-hand side of eq. (A14) 
we split the sum over % 2 and 23 into two pieces, one with % 2 < is and the second with i 2 > is- 
The first sum has i\ < % 2 < 23 and 14 < i 3 . It manifestly agrees with the right-hand side 
of eq. (A12), as given in eq. (A13), up to an overall factor of 1/2. The second sum has 
?4 < «3 < i 2 and i\ < i 2 . Relabelling the indices i\ <-> 24, i 2 <-> iz and reversing the order 
of the spinor string, we see that the second sum is precisely equal to the first sum. Adding 
the first and second sums together proves the identity (A12), which establishes the proper 
multi-particle factorization behavior (Al) of the amplitudes A s n (jjt) for p < n. 

Now consider the remaining cases where Kf_ n — > 0. These multi-particle poles are in fact 
the ones appearing in the I th term in the recursive construction (5.20). In eq. (5.28), the 
relevant l/s/... n poles are hidden in the terms with p = n — 1 in eq. (5.30) for S 2 - They can 
be found in the factor 

t2- { x-Y)ti- v \(P + l) + > = # 2 ...(i-i)#i...„ K) =~{ln) s,.. n , (A15) 

after using momentum conservation, if 2 ...(/-i) = —k\ — K[... n . Note also that in this limit, 
with P = —Ki... n , we have 

K (p+1) ... n = k n , (A16) 

(l-|^2...(i-i) = (1^>CP + I, (A17) 
Sl ... p = Si ... (n _D = (-P - k n f = (nP) [Pn] . (A18) 

Using these relations, the limiting behavior of the {l,p — n — 1) term in S 2 , as Kf„. n — > 0, is 

((I -I) I) 



S 2 (l,n-1) 



x 



(ln) 2 [nP] 2 (P (Z-l)) (PZ> 
((n — 1) n) 



x 



(IP) [Pn] (n(n-l))(ln) s,... n 
x(-l) x (lP) 2 [Pn] 2 (nl) 2 (jP) [Pn] (nl) 
IP) (P + \lF(l,n-l)] 2 \n-)(nl) 
(nP}[Pn] 

(IP) 3 UP) jiJ^m (nl) (n + \[F{l,n-l)f\p-) 

si... n ((l-l)P)(Pl)(nP)(Pl) [nP] 

(A19) 

If the second argument of jF(/,n — 1) in eq. (A19) were n instead of (n — 1), we would be 
done, as we could then use the same logic as in the case p < n treated earlier. But first we 
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have to show that 



(n+| \F(l,n-l)f\p-) = (n + \ n)] 2 \P~) . (A20) 

From the definition (A6), T{l,n) - F(l,n - 1) = hT^ZHj- The fir st T(l,n - 1) on 
the left-hand side of eq. (A20) can be replaced by jF(Z,n), because [nn\ = 0. The second 
one can be replaced as well, because the difference is proportional to Y^=i ( n ~\j = 
— (nP) [P P] = 0. This verifies the factorization behavior as Kf„, n — > 0. 
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